A0-A080  m  NAVAL  P0ST9AADUATE  SCHOOL  NONTEAEY  CA  f/%  fO/S 

DCS I ON  CONS IDEA AT IONS  OF  A  FLOW  STSTEN  FOA  TAANSIENT  LASEA  DISC— ETC 
AO#  TO  J  STAICKEA 

UNCLASSIFIED  NPS47-7*-###  hl 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS-  1963-A 


ADA  0  80  9  42 


NPS67- 79-008 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


Josef  Strieker 


August  1979 


Approved  for  public  release;  distribution  unlimited 

Prepared  for: 

Chief  of  Naval  Research 
Arlington,  VA  22217 


80  2  19  1 85 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


Rear  Admiral  T.  F.  Dedman  Jack  Borsting 

Superintendent  Provost 


The  work  reported  herein  was  supported  by  the  Foundation  Research 
Program  of  the  Naval  Postgraduate  School  with  funds  provided  by  the 
office  Chief  of  Naval  Research. 


Professor  Oscar  Biblarz  collaborated  closely  with  the  author  in  the 
preparation  of  this  report. 


Reproduction  of  all  or  part  of  this  report  is  authorized. 


This  report  was  prepared  by: 


Adjunct  Research  Professor 


of  Department  of  Aeronautics 


_ UNCLASSIFIED _ 

SfCUM'Tv  CL  ASSIFIC  ATtOM  Q>  Thi4  >»GI  (W*mm  Hhi  IaU»«<) _ 

ft  REPORT  DOCUMENTATION  PAGE  |  ■EFOMDco2Ke"sCVomi 

4  ujinr*  "| -  ✓  '  jj.  OOVT  ACCCSl'ON  NO. |  1  MCCIPICnT'S  CATALOG  NuaiC* 

NPS-*67- 79-008 


«.  TITLE  fan*  Mill)') 


jiPSIGN CONSIDERATIONS  OF  A^LOW  ^SYSTEM  FOR 

TRANS IE!nT  JASER  DISCHARGE  EXPERIMENTS  m, 

t _ ~  -*  •  "-s^ -  -  y 


i.  tvn*  or  wrorr  a  rcmoo  coveeeo 


1  Augue*  — - 15  Oc 


AuTmOACJ 


:ontmact  on  snant  nlnociw<) 


^JOSEF/STRICXER  / 

% _ oo»>onninu  iSHCTNiTATiOH  hmh  two  *oo»m 

NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  CA  93940 

n.  con t no l lino  orricc  nanc  and  aooncss 

NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  CA  93940 

14.  NONITONINO  AOCNCV  NADI  A  AOONCSVI*  MUarmtTl 


10.  NMOOIIAM  CLCNCNT.  NMOJCCT.  TAM 
AftCA  4  »«4K  UNIT  NUMIC4I 

61152N,RR  000-01-01 
N0001479WR90027 

™*  ’  12.  DUMf  OAT|  .  .  "" 

(//)  AuguBt 
u.  nubscw  orWAtfM 

31 _ 

i  Cmmtrmllln,  Othmm)  14.  tCCuni  .  T  CL  ASS.  (ml  Kim  r i*tn) 


Unclassified 


OCCL  AS4I  ric  ation/  DOWNOR  aoinc 
SCHEDULE 


10.  MSTNIOUTION  STATEMENT  (ml  Kim  Hmpmrl) 

Approved  for  Public  Release;  Distribution  unlimited 


0-z! 


I  17.  OtSTHiauTlOM  tTATCMCNT  (•<  <*•  Iff  »!••*  74,  It  *M—m$  ( 


(77)fth 


I  10.  lWeNLCMCNTANV  NOTH 

| Electrical  Discharges 
| Gasdynamics 
! Electrical  lasers 


Ludwieg  tubes 


j  12.  KIT  40*01  fCuNiN  mm  rmwmmmm  mlSm  II  iwmiimt  «N  lONIO  >r  Nw*  ‘ 


IO.  AOtTNACT  (Cmmllmmm  mm  rmmmtmm  ml  mm  ll'mmmmmmmr  mm*  immmlltf  Of  Mm*  mmmrnm) 

This  report  examines ’Resign  considerations^®^  a  flow  system  for  transient 
electrical  laser  experiments'.  After  a  brief  review  of  corona  discharge  gas  flow 
interaction  phenomena,  the  choice  of  the  size  and  operating  conditions  of  the 
discharge  cell  are  discussed.  Once  those  parameters  are  established,  the  flow 
system,  which  consists  essentially  of  a  Ludwieg  tube,  is  designed.  Attention 
is  paid  to  flow  starting  times,  to  the  available  time  of  steady  state  flow,  and  ta 
the  effects  of  the  boundary  layer  on  the  core  flow.  Suggested  experiments  and 
instrunentatidn  are  givenf  as-weliV  in  this  report. A _  I 


COITION  or  I  NOV  4t  12  OetOLCTC 
t/H  0  102*0 14*  440 1  I 


MCUMTT  CL  J 


UNCLASSIFIED 

40ICATION  or  TNI0  OAOO 


TABLE  OF  CONTENTS 


1.  Introduction  - 1 

2.  Test  Cell  Design  -  5 

3.  Ludwieg  Tube  - 9 

4.  Instrumentation  - 19 

5.  Figures  - 21 

6.  List  of  References  - 29 

7.  Distribution  List  -  31 


/ 

f 


Accession  For 


NT  IS  GRAAI 
DDC  TAB 
Unannounced 
Justification. 


Bv  - - - 

Distribution/ _ _ 

Avail nbllljT  Codes, 
Avail  and/or 
Dist  1  special 


R 


iii 


UST  OF  FIGURES 


FIGURE 

1  Schematic  of  Present  Discharge  Cell  and  Detail  of  Pin 

Region  -  21 

2  Transient  Flow  System  Schematic  -  22 

3  Suggested  Turbulence  Generating  Plates  Mounted  Close  to 

Electrode  Tips  — - 23 

4  Approximate  Wave  Diagram  for  the  Flow  in  a  Ludwieg  Tube. 

(a)  Upstream  Diaphragm  Location  -  24 

(b)  Downstream  Diaphragm  location  -  25 

5  Charge  Tube  Pressure,  P4»  as  a  function  of  Cell  Pressure, 

P2»  for  Different  Values  of  M2,  Y  =  1.4  - -  26 

6  Charge  Tube  Pressure,  P^,  as  a  function  of  Cell  Pressure, 

P2»  for  Different  Values  of  M2.  Y  =  1.67  -  27 

7  Vacuun  Tight  Compressor  - — - - -  28 


V 


abstract 


This  report  examines  design  considerations  of  a  flow  system  for 

t 

transient  electrical  laser  experiments.  After  a  brief  review  of  corona 
discharge  gas  flow  interaction  phenomena,  the  choice  of  the  size  and 
operating  conditions  of  the  discharge  cell  are  discussed.  Once  those 
parameters  are  established,  the  flow  system,  which  consists  essentially 
of  a  Ludwieg  tube,  is  designed.  Attention  is  paid  to  flow  starting 
times,  to  the  available  time  of  steady  state  flow,  and  to  the  effects 
of  the  boundary  layer  on  the  core  flow.  Suggested  experiments  and  in¬ 
strumentation  are  given,  as  well,  in  this  report. 
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I.  INTRODUCTION 


The  effect  of  gasdynamic  means  for  stabilizing  corona  discharges 
in  atmospheric  air  has  been  shown  to  be  significantly  important. 

By  choosing  the  right  turbulence  spectrum  and  flow  speed,  the  electri¬ 
cal  power  fed  into  a  diffuse  corona  discharge  has  been  improved  by  as 
much  as  500  times  above  the  maximum  no-flow  power.  The  influence  of 
turbulence,  which  is  more  substantial  at  the  higher  pressures  (-1  atm) , 
is  to  promote  better  mixing  and  homogeneity  of  the  gas.  Thus,  by  gas- 
dynamic  means,  molecular  gases  are  able  to  accept  higher  diffuse  dis¬ 
charge  power,  making  such  discharge  attractive  for  electric  gas  laser 
operation  because  of  the  promise  of  attractive  output  power  capabili¬ 
ties. 

The  temporal  non uniformities  of  the  gas  density  created  by  the 
turbulence  may  cause  phase  distortions  in  the  laser  output  beam,  and 
thus  deteriorate  the  laser  performance,  (due  to  the  effect  of  the  phase 
distortions  on  the  far  field  intensity  distribution) .  However,  turbu¬ 
lent  corona  discharge  studies  show  that  the  strongly  decaying  turbulence 
mainly  affects  the  discharge  at  the  region  close  to  anode  tips  ^ 
where  the  electric  field  is  strong  and  the  plasma  is  formed;  this  is 
the  corona  region  of  length  L^,  and  has  a  strong  nonuniform  electric 
field  Ea  (Fig.  1) .  The  external  zone  is  a  quasi-neutral  zone  which  has 
a  length  of  1^,  »  La,  and  an  electric  field  E^  which  is  considera¬ 

bly  lower  and  more  uniform  than  E  .  Thus,  by  exciting  the  turbulence 

A 

at  the  anode  tips  region  and  letting  the  turbulence  decay  and  become 
very  weak  at  the  external  zone,  where  laser  action  occurs,  the  density 
non uniformity  problem  can  be  reduced  to  minimum. 
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The  separation  of  the  discharge  region  into  two  main  regions  with  space 
scales  and  (where  «  L^)  and  electric  fields  E&  and  E^  (where 
Ea  >:>  Efa)  appears  to  be  very  useful  for  efficient  laser  operation. 

Those  regions  may  be  controlled,  almost  independently,  by  changing 
electrode  geometry  (anode  tips  radii,  electrode  separation,  electrode 
material,  etc.)  and  flow  characteristics.  It  is  possible  to  choose  the 
average  electric  field  in  the  external  region  so  that  the  E/N  value 
will  be  in  such  a  range  that  most  of  the  electron  energy  will  go  to 
excite  the  i?>per  laser  level;  as  for  example,  in  a  CO 2  laser  containing 
a  mixture  of  He/N2/co2  in  a  volumetric  ratio  of  3/2/1,  more  than  80%  of 
the  electron  energy  goes  into  the  CO^ (001)  and  N2 (v=l)  vibrational 
levels  for  the  value  of  E/N  s  2  x  lO-16  vcm2  .  In  this  respect, 

an  analogy  can  be  made  between  a  c.w.  electron-beam-sustained  electric 
discharge  and  a  c.w.  gas dynamic  stabilized  corona  discharge;  both  cases 
operate  in  the  Ionizer/Sustainer  mode  ^  ,  which  means  that  there  is  a 
separation  between  the  electric  field  that  produces  the  electrons  and 
the  sustained  discharge,  where  the  electrons  acquire  their  optimal 
temperature  Te  (Te  is  determined  by  E/N)  for  laser  action. 

In  this  report  we  propose  an  experimental  study  of  the  effects  of 
gas  dynamic  means  ( matin  ly  turbulence  and  flow  speed)  and  different  geo¬ 
metrical  parameters  of  the  electrodes  on  corona  discharges  for  laser 
applications . 

The  experimental  system  proposed  is  a  versatile  inexpensive  set-up, 
in  which  the  following  can  be  studied. 

Fluid  dynamics  measurements: 

Flow  rates 


Turbulence  spectrum 
Velocity  profiles 


Pressures 


Temperatures 


Electrical  measurements : 


V-I  characteristics 

probes  for  E,  Ne,  and  T  measurements 

e 

Optical  measurements: 

Photography 
Schlieren  photography 
In  terf e rograms 
Light  emission  spectroscopy 
Laser  performance  measurements : 

Small  signal  gain 
Output  power  measurements 
Parameters  to  be  studied; 

Turbulence  intensity  and  spectrum 
Flow  rates  (speed  and  pressure) 

Gas  mixtures 

Discharge  electrical  energy 
Distance  between  electrodes 
Electrodes  geometry 

Preionizing  devices,  like  u.v.  radiation 

To  generate  the  flow  through  the  corona  discharge ,  the  application 
(8) 

of  the  Ludwieg  tube  principle  has  been  considered  to  be  an  attractive 

proposition  because  of  its  simplicity,  low  cost  and  also  high  flow 

quality.  The  useful  flow  time  of  this  device  depends  mainly  on  the 

charge  tube  length  and  on  the  speed  of  sound  of  the  gas  in  the  tube  (see 

Section  III).  For  air  this  time  is  approximately  5  —  .  m3?c — —  .  For  a 

meter  length 

charge  tube  of  the  order  of  4m  in  length,  the  net  useful  flow  time  will 
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be  15-20  msec,  which  corresponds  to  -30  cavity  flow  times  (for  5  cm 

cavity  length  and  gas  speed  of  100  nv'sec) .  This  procedure  effectively 

simulates  cw  operation  without  requiring  the  heat  sink  electrodes, 

large  compressors ,  heat  exchanger  etc. ,  that  would  be  necessary  for 

true  c.w.  operation.  All  other  relevant  characteristic  times  are  much 

(2) 

shorter  than  the  flow  time 

Section  II  deals  with  the  basic  design  and  sizing  of  the  test  cell. 
Cnee  the  cell  dimensions,  the  pressure  of  the  gas  and  its  flow  speed  are 
determined,  the  Ludwieg  tube  flow  system  can  be  designed.  The  detailed 
considerations  of  the  design  are  given  in  Section  III . 
perimental  set-vqj  is  shown  schematically  in  Fig.  2. 


The  overall  ex- 
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II.  TEST  CELL  DESIGN 

The  test  section  design  is  basically  the  same  as  that  of  the  cell 

(2  3) 

used  previously  in  the  gasdynamic-discharge  interaction  studies  . 

Some  changes  have  to  be  made  due  to  the  use  of  this  cell  for  laser  ex¬ 
periments  and  as  being  a  part  of  a  different  flow  system. 

1.  Cell  Operating  Conditions 

Gas  mixtures  Different  gas  mixtures  are  going  to  be  used  in 
this  system:  N2,  He,  COj,  CO  ...  etc.  As  far  as  a  flow  system 
design  is  concerned,  the  difference  among  those  mixtures  is  their 
specific  heat  ratio  y  and  molecular  weight, 
will  consider  two  values  for  y : 
y  3  1.4  (diatomic  molecular  gas) 
y  =  1.67  (monatomic  gas) 

Gas  pressure  The  cell  will  be  operated  at  pressures  around  500 
torr  which  is  the  pressure  at  which  maximum  electrical  power  can  be 
dumped  into  the  corona  discharge  in  no  flowing  air  with  4cm  elec¬ 
trode  gap^ .  For  different  gases,  different  electrode  gaps  and  for 
flowing  gases  the  optimum  operating  pressure  varies,  thus  the  sys¬ 
tem  will  be  designed  to  operate  in  the  pressure  range  of  200  <  p  < 

800  Torr  (lower  and  higher  pressures  can  be  achieved  as  well) . 

Flow  speed  The  gasdynamics-discharge  interaction  becomes  effec¬ 
tive  at  flow  speeds  of  ~  100  m/sec ^  .  The  faster  the  flow  is,  the  more 
discharge  performance  improves.  The  cell  will  be  designed  to  have  a 
flow  of  about  Mach  number  M  =*  0.5. 
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2 .  Cell  dimensions 


Cross  section  dimensions: 

Cell  height:  d  =  2.22",  5.5  cm  . 

Cell  width:  b  =  6.5"  ,16  cm  , 

Dimension  along  the  flow  direction: 

Cell  length:  c  =  4"  ,10  cm  , 

Dimensions  d  and  c  are  the  same  as  those  in  Ref.  3;  however, 

the  cell  width,  b,  is  bigger  since  along  this  direction  the  small 

signal  gain  (s.s.g.)  is  measured  and  laser  cavity  will  be  aligned. 

(6  9) 

The  s.s.g.  reported  in  the  literature k  '  for  electrical  lasers, 
range  from  1%  to  4%  cm  16  cm  cavity  length  should  be  adequate  for 
gain  and  power  measurements. 

3.  Nozzle  throat 

The  cross  section  area  of  the  discharge  cell  is  S  =  b  x  d  = 

2 

88  cm  .  In  order  to  obtain  M  =  0.5  in  the  test  section,  the  cross 
section  area  A*,  at  the  throat  for  the  different  values  of  y,  are 
given  in  Table  I: 


y  =  1.4 

y  =  1.67 

A/A* 

1.34 

1.31 

A*  (cm2) 

66 

67 

h*(an) 

4.1 

4.2 

Table  I 


Table  I  give  also  the  values  of  throat  height  h*,  for  a  one 
dimensional  nozzle. 

Turbulence  generating  plates  have  to  be  inserted  ahead  of 
the  electrode  tips.  Since  the  effective  cross  section  of  those 
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plates  is  smaller  than  that  of  the  cell,  their  effective  cross 
section  should  exceed  the  values  of  A*  given  in  Table  X;  other¬ 
wise,  the  flow  will  choke  at  the  plate,  prior  to  reaching  the  test 
section. 

The  turbulent  plate  used  in  previous  experiments  had  an 
area  blockage  of  the  order  of  50%,  which  is  much  above  the  34% 
allowed  for  this  design.  A  new  design  and  location  of  those 
plates  have  to  be  considered.  One  way  of  improving  the  plate 
blockage  is  to  mount  the  turbulence  generators  closer  to  the 
electrode  tips.  This  will  enable  a  weaker  turbulence  generation 
,  (thus,  a  smaller  blockage)  since  the  decay  of  the  turbulence  at 
the  electrode  tips  will  be  insignificantly  small.  One  way  of 
doing  it  is  shown  in  Fig.  3. 

4.  Electrical  power  heat  addition 

The  effect  of  adding  heat  to  a  compressible  flow  is  to  send 
compression  waves  into  the  flow  which  change  flow  parameters;  heat 
addition  to  subsonic  flow  will  accelerate  the  flow  and  for  a  cri¬ 
tical  value  of  heat,  <3cr/  the  flow  will  be  choked. 

Q  ,  lc.r  , — i - &.+2»h?  -  x 

01  Cp  stag  2(y  +  1)M2  1  +  X-  M2 
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Typical  mass  flow  rates  through  the  cell: 


*  *  PAV  "  (r/5>t  Av 


for  N^: 


10 

2  x  88  x  2  x  10 

,8.3  x  107,  _ 

( - 55 - -U  300 


3  1  Kg/sec 


for  H  : 


ia  4 

i  2  x  88  x  2  x  10 

m  *  — - 7 -  3  0.14  Kg/sec 

,8.3  x  10  . 

( - : - )*■  300 


Typical  critical  heat  addition  values  are  given  in  Table  II: 


i 

r 


®cr 

qcr  k^kg/sec 

qtti  kw,For  typical 
m  values 

for  N2 

0.45 

140 

140 

for  B 

e 

0.4' 

620 

87 

The  power  supply  available  has  ~  1  kw  power  output  which  is 
very  low  power  compared  to  . 
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III.  LUDWIEG  TUBE 


Introduction 

The  Ludwieg  tube  flow  system  is  shown  schematically  in  Figure  2. 

The  Ludwieg  tube  consists  of  a  long  tube,  charge  tube,  where  the  gas 
is  initially  compressed  to  the  desired  charge  pressure,  A  tran¬ 

sition  section  channels  the  flow  from  the  circular  charge  tube  to 
the  rectangular  test  section  which  is  a  part  of  a  converging-diverging 
supersonic  nozzle.  The  gas  expanding  through  the  nozzle  is  discharged 
through  a  short  diffuser  into  a  dump  tank.  A  diaphragm  separates  the 
driver  and  driven  sections.  The  diaphragm  can  be  located  either  down¬ 
stream  or  upstream  of  the  nozzle. 

The  flow  is  initiated  by  the  rupture  of  the  diaphragm.  In  the 
case  of  an  upstream  diaphragm,  a  backwards  facing  centered  expansion 
fan  propagates  into  the  charge  tube.  This  fan  is  reflected  at  the 
closed  end  of  the  tube  and  returns  back  to  the  nozzle.  This  is  called 
the  first  flow  cycle;  a  simplified  wave  diagram  is  shown  in  Fig.  4(a) . 
Between  the  time  when  the  tail  of  the  expansion  wave  leaves  the  nozzle 
entrance  and  when  the  head  of  the  reflected  fan  reaches  the  nozzle 
again,  the  stagnation  conditions  ahead  of  the  nozzle  remain  constant  and 
a  clean  steady  flow  may  be  established  in  the  nozzle.  A  second  flow 
cycle  (and  others)  follows  the  first  one  but  the  stagnation  parameters 
are  different.  The  duration  of  the  steady  flow,  for  a  certain  gas,  is 
primarily  a  function  of  the  charge  tube  length  and  we  are  interested  to 
make  it  as  long  as  possible;  there  is  a  limit  on  the  length  of  the  tube 
set  by  the  growth  of  the  boundary  layer ^ 11 ' 12 ' .  The  boundary- layer 
growth  sets  also  a  limit  on  the  diameter  of  the  tube  and  on  the  flow  Mach 
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number,  M^,  in  the  charge  tube.  Subsequent  to  the  diaphragm  rupture, 
a  shock  wave  and  a  contact  surface  proceed  downstream  through  the 
driven  section.  Additional  waves  follow  the  contact  surface  to  adjust 
the  nozzle  test  section  pressure  to  the  freestream  pressure  correspond¬ 
ing  to  the  designed  Mach  number  expected  at  the  nozzle  exit  for  a  given 
(14) 

nozzle  .  The  steady  supersonic  flow  in  the  nozzle  is  not  established 
until  all  those  waves  are  swept  out  the  nozzle.  The  time  required  to 
get  the  flow  started,  that  is,  the  time  required- to  achieve  a  steady 
flow  from  rest,  in  the  supersonic  nozzle,  is  a  very  important  performance 
criterion.  It  is  important  to  minimize  the  start  time  in  order  to  maxi¬ 


mize  the  duration  of  the  steady  state  conditions  for  a  fixed  charge  tube 
length(14'15). 


For  downstream  diaphragm  location,  the  wave  diagram  is  different 
from  the  case  where  the  diaphragm  is  located  upstream,  a  simplified  x  -  t 
diagram  is  shown  in  Fig.  4(b) .  As  soon  as  the  local  sonic  velocity  is 
reached  at  the  nozzle  throat,  no  further  expansion  waves  can  pass  through 
into  the  tube  and  an  expansion  fan  of  finite  width  passes  down  the  tube. 

For  an  upstream  diaphragm  location,  the  starting  times  at  the 
nozzle  exit  are  reported  to  be  shorter  than  those  obtained  with  the 
downstream  diaphragm  location .  However,  the  disadvantage  of  the  up¬ 
stream  location  diaphragm  is  that  when  breaking  the  diaphragm .particles 
might  tear  off  and  disturb  the  flow  in  the  test  section. 

2.  Charge  tube  diameter  and  length 

Economic  considerations  dictate  that  a  minimum  charge  tube  diameter 
be  used  for  a  given  test  section  size.  This  results  in  a  relatively 
small  contraction  (small  ratio  between  nozzle  throat  area  to  charge  tube 


( 


area)  which  requires  high  Mach  number,  M3,  at  the  charge  tube.  This 

will  cause  the  boundary  layer  to  have  a  dominant  effect  on  nozzle  flow. 

It  is  suggested  to  have  the  I . 0.  of  the  charge  tube  to  be  the  same 

as  the  dimension  b  of  the  test  cell,  i.e.  16  cm.  In  this  case  A./A* 

4 

and  are  given  in  Table  III. 


Y  =  1.4 

Y  =  1.67 

a4/a* 

3.05 

3.00 

M3 

0.2 

0.19 

Table  III 


If  a  tube  with  a  larger  diameter  is  available  (like  20cm) ,  it  will  im¬ 
prove  the  performance  of  the  system,  since  M3  will  be  smaller.  It  is 
not  recommended  to  use  a  plenum  chamber  upstream  of  the  nozzle  inlet 
(to  slow  and  mix  the  core  and  boundary  layer  flows  emerging  from  the 
pipe)  since  it  will  be  shown  that  for  M3  *  0.2  the  disturbances  intro¬ 
duced  by  the  boundary  layer  growth  are  small  while  the  starting  time  is 
going  to  be  much  longer  by  introducing  the  plenvm  chamber. 

6*/d  is  given  as  a  function  of  instananeous  aspect  ratio  of  the 

(P  d)1/5 

flow  s/d,  P.d  and  M  ,  in  the  form  of  5*/d  -  by  Russel ,  based 

(s/d) 4/5 

on  Becker's  model  for  5*/d^^  ,  where  6*  is  the  bound  jury  layer  displace¬ 
ment  thickness,  at  a  distance  s  from  the  expansion  wave  head,  and  d  is 
the  charge  tube  diameter. 


5*  <V)V5 

for  M.  -  0.2  -  - 

3  d  (s/d) 4/5 


5  x  10  4  (atm  -  m)1^ 


At  the  end  of  the  test  time,  s  “  2L  and  the  thickness  of  the  boun¬ 
dary  layer  is  5  *  106*.  For  P4  -  1  atm,  L  -  3  m  and  d  -  0.16  m. 
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we  get  6*/d  =*  1.3%,  5  “  0.02  m  and  6/d  *13%. 

The  perturbations  in  the  core  flow  due  to  the  boundary  layer, 
when  the  test  period  is  over,  are  (for  the  same  parameters  used  for 

6* 

— —  calculations)  : 

<z 


-AP3 

(P.d) 1/5 

4  .  i  3 

y  1Q-4 

P3 

*  4*  J 

(L/d) 4/5 

X  XU 

P3 

(P4d)1/S 

x  10~3 

*  fa. 

V 

3 

(L/d)4/5  3*2  ' 

v3 

at3 

Y  -  1  AP3 

^i.X6 

rr x  6 

.5  x  10~3- 

T3 

Y  P3 

6.5  x  10 


-3 


4.8  x  10 


-2 


which  are  relatively  very  small  disturbances . 

The  effects  of  the  above  described  core  perturbations  on  the  down¬ 
stream  flow  were  discussed  by  Russell  .  He  evaluated  the  boundary 
layer  along  the  walls  of  a  set  of  specified  nozzles.  The  results  for 
different  values  of  5*/d,  and  M  (Mach  number  in  the  nozzle)  are 
presented  in  Fig.  7  of  his  paper;  he  presents  the  effective  area  - 

as  a  function  of  6*/d  for  different  values  of  M,  and 
A/A*  3 

Mf where  (A/A*)#f£  is  the  test  section  -  to  -  throat  area  ratio  with 
boundaxy-layer  correction. 

A/A*  is  the  geometric  area  ratio,  with  no  corrections. 
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For  the  value  of  P,d  *  0,5  atm-m,  Russell  found  that  with 
4 

10~4  <  <5*/d  <  the  effective  area  change  parameter  has  a  finite 

value  determined  by  the  normal  steady  boundary  layer  correction  of 

the  nozzle  itself  (in  our  case  P  d  B  2  x  0.16  =  0.3.  There  is  a 

o 

weak  dependence  on  pressure  level  (P^d)  )  • 

TO  summarize : 

for  M  =■  0.5  Mach  nunfcer  flow  in  the  test  section/  d«16cm,  L=3m 
and  m3  *  0.2 

<5*/d  -  1.3%  t  6/d  ■  13% 

AP 

—  -  0.6%  ,  -  « 

(A/A*)eff  -  A/A*  ^  ^  (contribution  from  charge  tube 
A/A*  boundary  layer) 


3.  Test  time 

The  test  time,  (the  Him  taken  for  the  head  of  the  reflected 
expansion  wave  to  reach  back  the  throat)  for  small  values  of  M3  is 
given  approximately  by  2L/a^,  where  a4  is  the  sound  speed  of  the 
undisturbed  gas  in  the  charge  tube. 

The  time  available  falls  only  slightly  with  increasing  the  Mach 

.  „  (17) 

number  M3 


t/t  ■ 
o 

where  t  ”  L/a4«  For 
t  and  t  are  given  in 


■T  +  1 

2  {1  ,  y  -  i  M  }2(y  -  i) 

1  +  M,  U  +  **3^ 


M3  «  0.2  and  L  =  3m,  the  calculated  values  of 


Table  IV. 
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Y  *  1.4 

Y  -  1.67 

a^ (m/sec) 

350 

1000 

t  (msec) 
o 

8.55 

3 

t/fco 

1.88 

1.9 

t(msec) 

16 

5.7 

Table  XV 


4.  Starting  times 

The  starting  time  of  supersonic  flow  at  a  given  location  in  the 
flow  system  is  defined  as  the  time  elapsed  between  the  arrival  of  the 
head  of  the  expansion  fan  (downstream  diaphragm) ,  or  shock  wave 
(upstream  diaphragm)  and  the  instance  where  the  pressure  becomes  time  in' 
dependent,  i.e.,  when  steady  flow  is  established.  For  an  upstream 
diaphragm  location,  the  starting  times  at  the  nozzle  exit  are  reported 
to  be  shorter  than  those  obtained  with  the  downstream  diaphragm  loca- 


a.  Downstream  diaphragm  location: 

According  to  ref.  15,  (they  used  a  Ludwieg  Tube  with  a  charge  tube 
tube  of  67C  cm  in  length  and  X.D.  of  13.4  cm)  the  starting  times  at  the 
nozzle  exit  are  independent  of  the  initial  downstream  pressure  p  , 
except  in  the  range  above  performance  limit.  Starting  times  increase 
with  diffuser  length.  Die  starting  times  upstream  of  the  nozzle  throat 
do  not  depend  noticeably  on  diffuser  length,  and  are  of  the  order  of 
0.5  msec,  where  the  starting  times  in  the  diverging  part  of  the  nozzle 
are  of  the  order  of  6  msec.  The  experimental  results  are  presented  in 


dimensionless  starting  time  ts/T£  as  a  function  of  dimensionless  axial 
distance  x/D ,  where  is  defined  as  £  *  </R*h*,  R*  is  the  radius 

of  curvature  at  the  nozzle  throat,  h*  the  nozzle  height  and  a*  the 
speed  of  sound  at  the  throat.  At  the  subsonic  nozzle  section  tg/x^  s  3 
at  the  nozzle  exit  t  /t_  a  40.  In  our  case,  the  starting  time  in  the 

S  £ 

subsonic  section  will  be  tg3  1  msec  and  in  the  supersonic  section 

t  3  8  msec,  for  N_.  For  other  gases,  the  starting  time  scales  in- 
s  2 

versely  to  a*. 

b.  Upstream  diaphragm  location 

The  starting  times,  at  the  exit  of  the  nozzle,  increase  with 
initial  downstream  pressure,  P^,  and  increase  also  with  diffuser  length. 
The  starting  times  are  of  the  order  of  2  msec  at  nozzle  exit  and  P^/P^ 
should  be  at  least  of  the  order  of  10  in  order  to  start  the  flow^*  . 
Summary: 

It  is  recommended  to  use  the  downstream  diaphragm  location, 

mainly  because  the  starting  time  in  the  subsonic  portion  of  the  nozzle 

is  very  short  ( -  1  msec)  and  the  troubles  of  having  diaphragm  particles 

in  the  flow  are  eliminated.  For  supersonic  experiments,  it  might  be  an 

advantage  to  have  the  diaphragm  located  upstream.  The  diaphragm  holder 

unit  cam  be  designed  so  that  there  will  be  a  possibility  to  mount  it 

either  downstream  or  upstream. 

4.  Charge  tube  pressures  and  temperatures 

The  pressures  and  temperatures,  P  and  T  ,  in  the  charge  tube,  prior 

4-4 

to  diaphragm  rupture,  am  a  function  of  pressure  ?2  and  temperature 
respectively,  where  the  Mach  number  of  the  flow  is  M^,  is  given  by : 
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4  st  4  3  st  (1  +  - — M  )  1  + 


P_  T_  - *-rr 

-i-  t-2*y  - 1 


2  1 

For  y  *  1.4 ,  M-  *  0.2  -*-=■  =  0.93  - - - - 

4  1  +  M2 


2  1 

For  Y  *  1.67,  M,  =  0.2  -*•-=-  =  0.95  - = — - 

3  T,  .  Y  -  1  2 

4  1  +  2~  M2 

The  corresponding  values  of  P4/P2  can  be  calculated  from  Eq.  (2) . 

as  a  function  of  P2  for  different  values  of  the  is  she 
in  Fig.  5  and  Fig.  6  for  y  ■  1.4  and  y  =  1.67  respectively.  The 
values  of  T2,  for  Y  -  1.4  and  y  =  1.67  and  different  values  of  M2  a 
given  in  Table  V  assuming  T4  =  300°k : 


0 

0.5 

1 

1.5 

2 

280 

266 

233 

193 

155 

285 

263 

213 

162 

122 

Table  V 


The  pressure,  P^,  in  the  part  of  the  system  downstream  the  diaphragm 
should  be  of  the  order  of  P4/100  (to  get  best  start  times  without  a 
diffuser) ;  i.e.  P,  *  1  torr  (see  Sec.  5) . 


The  starting  times,  as  discussed  before,  become  longer  with  longer 
diffusers,  both  for  downstream  and  upstream  diaphragm  locations. 


For  a  fixed  charge  tube  pressure/  and  back  pressure  being  in¬ 
creased  gradually,  the  flow  at  the  diffuser  end  ceases  to  start  for  a 
certain  value  of  P^  while  further  increasing  P^,  the  supersonic  flow  at 
the  nozzle  exit  barely  starts.  Beyond  these  critical  conditions,  the 
flow  at  the  nozzle  exit  becomes  unstable.  The  limiting  experimental 
downstream  pressure  radios,  P^/P^  were  measured  by  Merkli  and  flbuaf , 
both  for  upstream  and  downstream  diaphragm  locations,  as  a  function  of 
diffuser  length  L/D. 

It  was  found  that  (P^/P^ cr  for  diffuser  starting  is  -0.11  inde¬ 
pendent  on  L/D. 

(P1/P4) cr  for  nozzle  supersonic  flow  increases  with  L/D  and  reaches 

a  maximum  at  about  L/D  *  6  where  (P,/P.)  *  0.22. 

1  4  cr 

That  means  that  within  the  available  time,  the  pressure  in  the  dump 
tank  should  not  increase  above  0.1  x  P^;  this  limiting  pressure  defines 
the  volume  of  the  dump  tank. 

The  rate  of  pressure  rise  inside  the  tank:  The  continuity  dictates 

apd 

that  in  the  dun?)  tank  +  m^  *0  where  is  the  tank's  volume,  p^ 

is  the  density  and  m^  is  the  total  mass  flow  into  the  tank.  By  differen¬ 
tiating  the  equation  of  state,  P^  =  pd  R  T^,  with  respect  to  time  and 
using  the  continuity  equation,  we  get 


RT 


d  at 


■d  ■  "a  ■  p3  A3  v3  * 


P3  A3M3a3 
RT- 


AP, 


TdP3A3M3a3 

T3Vd 


17 


Where  t  is  the  available  testing  time.  Since  initially  P  =  p 

d  J 

0.1P3  we  may  write  AP^  =  p^j  from  the  last  equation: 

Td  P3 

=  (p”)A3M3*BT  -  2  x  10  x  nx82x0.2  x  33000  x  20  x  10-3 

3  d 

V.  *  5.3  x  105  on3 
a 


It  is  suggested  to  install  a  valve  \q>stream  the  diapragm,  as 
shown  in  Fig.  3.  After  each  run,  before  the  ruptured  diaphragm  is 
taken  off,  this  valve  will  be  closed  and  the  gas  in  the  system  will 
be  saved. 

A  more  efficient  way  of  saving  gas  will  be  to  use  a  compressor 
to  compress  the  gas  in  the  dump  tank  and  recirculate  it  back  to  the 
charge  tube.  The  compressor  should  be  vacuum  tight.  A  simple  and 
inexpensive  compressor  is  suggested  in  Fig.  7. 


<< 
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IV.  INSTRUMENTATION 


Pressure  gauges 

0-5  atm.  range  for  measuring  P4 
0-100  torr  - 

0-760  torr  Pressure  gauges  for  gas  mixtures  preparation  . 

0-2  atm 

0-100  torr  for  measuring  P1 

For  measuring  and  P3  during  the  experiment  itself,  after 
diaphragm  is  broken,  we  need  piezoelectric  pressure  transducers  with 
time  response  of  the  order  of  tens  microseconds,  in  the  range  of 
0-2  atm. 


Velocity  profiler  and  turbulence  spectrum 

For  transient  flow  speed  and  turbulence  measurements,  a  laser- 
doppler  anemometer  or  a  hot-wire  anemometer  can  be  used.  The  avail¬ 
able  steady  flow  time,  10-20  msec,  might  be  too  short  for  appropriate 
averaging  the  turbulence  spectrum,  mainly  in  the  low  frequency  regime. 
Therefore  it  will  be  necessary  to  run  the  flow  in  a  continuous  mode  (with 
compressed  air)  for  turbulence  spectrum  measurements.  The  transient  flow 
fluctuations  may  not  be  averageable  but  their  ability  to  stabilize  the 
discharge  may  surpass  that  of  steady  flow  in  which  case  other  measurement 
forms  should  be  sought. 

Voltage  current  measurements;  probes 

The  voltage  and  current  have  to  be  recorded  during  the  experiment. 
The  simplest  way  of  doing  it  would  be  by  using  an  oscilloscope. 

The  probes  used  for  measuring  plasma  parameters  must  be  provided 
with  automatic  bias  sweep,  so  that  the  characteristics  can  be 
accomplished  within  the  flow  time. 
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Lasing  diagnostics 


A  probe  laser  for  s.s.g.  measurements. 

A  fast  response  detector  (microseconds  range) . 


Power  meter. 


Fig  2.  TRANSIENT  FLOW  SYSTEM  SCHEMATIC. 


TO  ELECTRODE  TIP. 


EXPANSION  FAN 


(1) 


L  SHOCK  WAVE 
L  CONTRACT  SURFACE 

UPSTREAM  FACING  WAVE 


Fig.  4(a)  APPROXIMATE  WAVE  DIAGRAM  FOR  THE  FLOW  IN  A 
LUDWIEG  TUBE.  —  UPSTREAM  DIAPHRAGM  LOCA- 


EXPANSION  FAN  \  \ 

\  SHOCK  WAVE 
CONTACT  SURFACE 


Fig.  4(b)  APPROXIMATE  WA/E  DIAGRAM  FOR  THE  FLOW  IN  A 

LUOWIEG  TUBE.- DOWNSTREAM  OIAPHRAGM  LOCATION. 


P2  (torr) 


Fig.  5  CHARGE  TUBE  PRESSURE,  P,  ,  AS  A  FUNCTION  OF  CELL 

4 

PRESSURE,  P2,  FOR  DIFFERENT  VALUES  OF  Mg.  y  «1.4 


P,  (torr) 
4 


1800 


Y  =1.67 


M2=  2 


M2*l.5 


1600 


M_  *  1 
2  > 


1400 


1200 


1000 


M  *  0.5 
2 


D  100  200  300  400  500  600 

P2  (torr) 

Fig.  6  CHARGE  TUBE  PRESSURE,  ,  AS  A  FUNCTION  OF  CELL  PRESSURE, 
FOR  DIFFERENT  VALUES  OF  y  *1.67 
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